This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.3793103

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.05 while 0.95 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_black

11

18

2

0.879

0.917

0.558

0.988

1.011

0.113

-0.337

CSEO

p_yellow

13

1

1

0.420

0.363

0.126

0.980

1.057

0.225

0.254

NSEI

p_yellow

13

1

2

0.424

0.371

0.133

0.978

1.054

0.223

0.238

NSEI

p_yellow

13

10

1

0.423

0.369

0.126

0.980

1.051

0.225

0.240

NSEI

p_yellow

13

10

2

0.424

0.370

0.126

0.978

1.062

0.225

0.241

NSEI

p_yellow

13

11

1

0.416

0.362

0.124

0.979

1.053

0.227

0.237

NSEI

p_yellow

13

11

2

0.420

0.361

0.123

0.978

1.051

0.225

0.262

NSEI

p_yellow

13

12

1

0.422

0.371

0.126

0.976

1.055

0.223

0.228

NSEI

p_yellow

13

12

2

0.426

0.380

0.127

0.977

1.052

0.224

0.206

NSEI

p_yellow

13

13

1

0.422

0.371

0.131

0.977

1.055

0.222

0.230

NSEI

p_yellow

13

13

2

0.427

0.374

0.128

0.981

1.049

0.222

0.238

NSEI

p_yellow

13

14

1

0.418

0.367

0.117

0.977

1.056

0.225

0.227

NSEI

p_yellow

13

14

2

0.427

0.372

0.131

0.975

1.057

0.223

0.247

NSEI

p_yellow

13

15

1

0.416

0.362

0.128

0.978

1.060

0.223

0.242

NSEI

p_yellow

13

15

2

0.419

0.359

0.128

0.982

1.055

0.223

0.269

NSEI

p_yellow

13

16

1

0.421

0.366

0.125

0.979

1.052

0.224

0.245

NSEI

p_yellow

13

16

2

0.422

0.370

0.124

0.976

1.054

0.223

0.233

NSEI

p_yellow

13

17

1

0.417

0.361

0.125

0.976

1.054

0.223

0.251

NSEI

p_yellow

13

17

2

0.428

0.377

0.125

0.977

1.053

0.223

0.229

NSEI

p_yellow

13

18

1

0.421

0.368

0.119

0.979

1.056

0.225

0.236

NSEI

p_yellow

13

18

2

0.421

0.371

0.119

0.978

1.050

0.224

0.223

NSEI

p_yellow

13

19

1

0.420

0.367

0.117

0.975

1.050

0.224

0.237

NSEI

p_yellow

13

19

2

0.426

0.374

0.132

0.979

1.054

0.223

0.233

NSEI

p_yellow

13

2

1

0.423

0.370

0.121

0.978

1.057

0.225

0.236

NSEI

p_yellow

13

2

2

0.424

0.370

0.127

0.980

1.060

0.224

0.241

NSEI

p_yellow

13

20

1

0.418

0.363

0.124

0.980

1.052

0.224

0.246

NSEI

p_yellow

13

20

2

0.423

0.365

0.130

0.979

1.055

0.223

0.260

NSEI

p_yellow

13

21

1

0.419

0.364

0.126

0.977

1.057

0.223

0.247

NSEI

p_yellow

13

21

2

0.423

0.365

0.133

0.980

1.052

0.224

0.259

NSEI

p_yellow

13

22

1

0.418

0.362

0.121

0.979

1.057

0.223

0.251

NSEI

p_yellow

13

22

2

0.420

0.367

0.124

0.976

1.053

0.222

0.239

NSEI

p_yellow

13

23

1

0.388

0.333

0.100

0.972

1.060

0.228

0.241

NSEI

p_yellow

13

23

2

0.472

0.429

0.155

0.983

1.040

0.217

0.198

NSEI

p_yellow

13

24

1

0.387

0.328

0.100

0.974

1.063

0.227

0.260

NSEI

p_yellow

13

24

2

0.470

0.426

0.155

0.982

1.052

0.218

0.202

NSEI

p_yellow

13

25

1

0.386

0.332

0.100

0.975

1.058

0.224

0.241

NSEI

p_yellow

13

25

2

0.470

0.431

0.152

0.982

1.051

0.218

0.179

NSEI

p_yellow

13

26

1

0.382

0.320

0.103

0.971

1.053

0.227

0.274

NSEI

p_yellow

13

26

2

0.466

0.428

0.152

0.981

1.050

0.217

0.175

NSEI

p_yellow

13

27

1

0.369

0.309

0.100

0.969

1.058

0.221

0.272

NSEI

p_yellow

13

27

2

0.462

0.419

0.154

0.981

1.042

0.215

0.200

NSEI

p_yellow

13

28

1

0.383

0.331

0.100

0.960

1.043

0.215

0.242

NSEI

p_yellow

13

28

2

0.465

0.429

0.155

0.975

1.038

0.210

0.171

NSEI

p_yellow

13

29

1

0.362

0.314

0.103

0.926

1.036

0.202

0.238

NSEI

p_yellow

13

29

2

0.456

0.423

0.151

0.953

1.035

0.201

0.164

NSEI

p_yellow

13

3

1

0.423

0.371

0.117

0.977

1.061

0.224

0.232

NSEI

p_yellow

13

3

2

0.423

0.369

0.126

0.980

1.054

0.224

0.242

NSEI

p_yellow

13

4

1

0.419

0.365

0.119

0.979

1.055

0.225

0.240

NSEI

p_yellow

13

4

2

0.424

0.368

0.123

0.980

1.053

0.226

0.248

NSEI

p_yellow

13

5

1

0.420

0.367

0.124

0.976

1.051

0.225

0.235

NSEI

p_yellow

13

5

2

0.422

0.367

0.127

0.978

1.043

0.226

0.244

NSEI

p_yellow

13

6

1

0.417

0.361

0.119

0.982

1.047

0.226

0.248

NSEI

p_yellow

13

6

2

0.425

0.375

0.120

0.979

1.054

0.225

0.222

NSEI

p_yellow

13

7

1

0.415

0.361

0.124

0.979

1.058

0.226

0.239

NSEI

p_yellow

13

7

2

0.426

0.373

0.125

0.979

1.063

0.221

0.239

NSEI

p_yellow

13

8

1

0.423

0.371

0.122

0.979

1.053

0.223

0.233

NSEI

p_yellow

13

8

2

0.424

0.371

0.125

0.978

1.053

0.224

0.236

NSEI

p_yellow

13

9

1

0.419

0.367

0.124

0.980

1.053

0.222

0.234

NSEI

p_yellow

13

9

2

0.421

0.366

0.128

0.977

1.051

0.224

0.245

NSEI

p_yellow

14

1

1

0.579

0.559

0.249

0.984

1.018

0.190

0.106

NSEO

p_yellow

14

1

2

0.618

0.610

0.292

0.989

1.017

0.181

0.044

NSEO

p_yellow

14

10

1

0.579

0.563

0.262

0.983

1.012

0.186

0.086

NSEO

p_yellow

14

10

2

0.617

0.606

0.296

0.987

1.012

0.181

0.061

NSEO

p_yellow

14

11

1

0.578

0.563

0.251

0.983

1.013

0.187

0.080

NSEO

p_yellow

14

11

2

0.621

0.616

0.275

0.989

1.010

0.184

0.027

NSEO

p_yellow

14

12

1

0.581

0.569

0.254

0.984

1.012

0.187

0.064

NSEO

p_yellow

14

12

2

0.621

0.614

0.288

0.989

1.012

0.180

0.039

NSEO

p_yellow

14

13

1

0.580

0.569

0.244

0.983

1.014

0.189

0.058

NSEO

p_yellow

14

14

1

0.580

0.571

0.250

0.986

1.018

0.186

0.048

NSEO

p_yellow

14

14

2

0.619

0.609

0.283

0.987

1.010

0.181

0.055

NSEO

p_yellow

14

15

1

0.578

0.561

0.252

0.985

1.017

0.189

0.090

NSEO

p_yellow

14

15

2

0.620

0.615

0.289

0.988

1.010

0.181

0.028

NSEO

p_yellow

14

16

1

0.580

0.562

0.253

0.984

1.010

0.190

0.095

NSEO

p_yellow

14

17

1

0.577

0.563

0.249

0.984

1.013

0.190

0.074

NSEO

p_yellow

14

17

2

0.615

0.607

0.275

0.987

1.015

0.186

0.043

NSEO

p_yellow

14

18

1

0.573

0.554

0.250

0.984

1.013

0.188

0.101

NSEO

p_yellow

14

18

2

0.618

0.611

0.288

0.989

1.015

0.180

0.039

NSEO

p_yellow

14

19

1

0.580

0.560

0.241

0.987

1.017

0.188

0.106

NSEO

p_yellow

14

19

2

0.622

0.613

0.294

0.988

1.012

0.181

0.050

NSEO

p_yellow

14

2

1

0.582

0.569

0.254

0.986

1.013

0.188

0.069

NSEO

p_yellow

14

20

1

0.579

0.564

0.258

0.986

1.012

0.187

0.080

NSEO

p_yellow

14

21

1

0.578

0.565

0.255

0.986

1.016

0.188

0.069

NSEO

p_yellow

14

21

2

0.617

0.609

0.283

0.989

1.014

0.183

0.044

NSEO

p_yellow

14

22

1

0.583

0.569

0.259

0.985

1.016

0.185

0.076

NSEO

p_yellow

14

22

2

0.620

0.609

0.298

0.988

1.010

0.181

0.061

NSEO

p_yellow

14

23

2

0.622

0.612

0.287

0.990

1.011

0.180

0.055

NSEO

p_yellow

14

24

1

0.608

0.599

0.292

0.983

1.011

0.174

0.052

NSEO

p_yellow

14

24

2

0.621

0.613

0.293

0.988

1.014

0.183

0.044

NSEO

p_yellow

14

25

1

0.594

0.579

0.272

0.986

1.016

0.179

0.084

NSEO

p_yellow

14

25

2

0.616

0.606

0.280

0.989

1.012

0.185

0.054

NSEO

p_yellow

14

26

1

0.588

0.574

0.260

0.987

1.013

0.186

0.075

NSEO

p_yellow

14

26

2

0.618

0.612

0.280

0.990

1.014

0.185

0.032

NSEO

p_yellow

14

27

1

0.602

0.593

0.264

0.986

1.012

0.182

0.049

NSEO

p_yellow

14

27

2

0.622

0.614

0.277

0.989

1.015

0.181

0.044

NSEO

p_yellow

14

28

1

0.584

0.571

0.260

0.983

1.013

0.181

0.072

NSEO

p_yellow

14

28

2

0.615

0.608

0.269

0.987

1.010

0.184

0.038

NSEO

p_yellow

14

29

1

0.584

0.576

0.242

0.971

1.011

0.182

0.044

NSEO

p_yellow

14

29

2

0.614

0.609

0.274

0.978

1.011

0.179

0.028

NSEO

p_yellow

14

3

1

0.577

0.565

0.248

0.986

1.016

0.192

0.063

NSEO

p_yellow

14

3

2

0.625

0.619

0.292

0.989

1.012

0.178

0.034

NSEO

p_yellow

14

4

1

0.579

0.566

0.243

0.984

1.012

0.190

0.069

NSEO

p_yellow

14

4

2

0.620

0.612

0.284

0.988

1.014

0.182

0.044

NSEO

p_black

14

4

1

0.920

0.947

0.682

0.993

1.013

0.083

-0.324

NSEO

p_yellow

14

5

1

0.579

0.568

0.250

0.985

1.013

0.191

0.058

NSEO

p_yellow

14

5

2

0.620

0.611

0.296

0.988

1.015

0.180

0.050

NSEO

p_yellow

14

6

1

0.577

0.561

0.254

0.984

1.011

0.187

0.086

NSEO

p_yellow

14

6

2

0.623

0.620

0.291

0.988

1.011

0.182

0.017

NSEO

p_yellow

14

7

1

0.576

0.560

0.248

0.986

1.014

0.190

0.084

NSEO

p_yellow

14

7

2

0.619

0.613

0.287

0.987

1.015

0.184

0.033

NSEO

p_yellow

14

8

2

0.621

0.614

0.283

0.990

1.013

0.182

0.038

NSEO

p_yellow

14

9

1

0.578

0.560

0.256

0.985

1.014

0.188

0.096

NSEO

p_yellow

14

9

2

0.620

0.609

0.295

0.988

1.015

0.181

0.061

NSEO

p_pelagic

16

1

1

0.370

0.363

0.120

0.656

1.052

0.145

0.048

SSEO

p_pelagic

16

1

2

0.349

0.340

0.108

0.638

1.052

0.142

0.063

SSEO

p_pelagic

16

10

1

0.372

0.364

0.130

0.664

1.036

0.141

0.057

SSEO

p_pelagic

16

10

2

0.353

0.344

0.117

0.641

1.048

0.140

0.064

SSEO

p_pelagic

16

11

1

0.375

0.369

0.132

0.670

1.045

0.143

0.042

SSEO

p_pelagic

16

11

2

0.357

0.349

0.115

0.643

1.054

0.140

0.057

SSEO

p_pelagic

16

12

1

0.377

0.371

0.134

0.661

1.039

0.139

0.043

SSEO

p_pelagic

16

12

2

0.356

0.348

0.114

0.630

1.053

0.136

0.059

SSEO

p_pelagic

16

13

1

0.380

0.373

0.141

0.663

1.039

0.139

0.051

SSEO

p_pelagic

16

13

2

0.361

0.351

0.127

0.641

1.051

0.135

0.074

SSEO

p_pelagic

16

14

1

0.365

0.358

0.137

0.638

1.039

0.132

0.053

SSEO

p_pelagic

16

14

2

0.358

0.350

0.131

0.626

1.039

0.132

0.060

SSEO

p_pelagic

16

15

1

0.373

0.364

0.141

0.654

1.028

0.133

0.068

SSEO

p_pelagic

16

15

2

0.359

0.351

0.128

0.640

1.051

0.132

0.061

SSEO

p_pelagic

16

16

1

0.369

0.357

0.137

0.649

1.035

0.136

0.088

SSEO

p_pelagic

16

16

2

0.362

0.347

0.140

0.638

1.026

0.132

0.114

SSEO

p_pelagic

16

17

1

0.421

0.417

0.188

0.683

1.029

0.130

0.031

SSEO

p_pelagic

16

17

2

0.378

0.368

0.145

0.654

1.036

0.134

0.075

SSEO

p_pelagic

16

18

1

0.406

0.400

0.172

0.672

1.032

0.133

0.045

SSEO

p_pelagic

16

18

2

0.374

0.369

0.149

0.651

1.027

0.133

0.038

SSEO

p_pelagic

16

19

1

0.400

0.393

0.175

0.662

1.024

0.128

0.055

SSEO

p_pelagic

16

19

2

0.379

0.371

0.160

0.649

1.022

0.129

0.062

SSEO

p_pelagic

16

2

1

0.366

0.361

0.119

0.658

1.051

0.143

0.035

SSEO

p_pelagic

16

2

2

0.348

0.341

0.107

0.633

1.056

0.141

0.050

SSEO

p_pelagic

16

20

1

0.400

0.395

0.178

0.661

1.014

0.127

0.039

SSEO

p_pelagic

16

20

2

0.387

0.380

0.161

0.647

1.018

0.128

0.055

SSEO

p_pelagic

16

21

1

0.388

0.380

0.168

0.659

1.019

0.127

0.063

SSEO

p_pelagic

16

21

2

0.380

0.373

0.161

0.645

1.022

0.126

0.056

SSEO

p_pelagic

16

22

2

0.387

0.383

0.167

0.645

1.011

0.123

0.033

SSEO

p_pelagic

16

3

1

0.368

0.361

0.122

0.659

1.048

0.143

0.049

SSEO

p_pelagic

16

3

2

0.348

0.338

0.110

0.627

1.064

0.141

0.071

SSEO

p_dsr

16

31

1

0.402

0.401

0.376

0.433

1.665

0.015

0.066

SSEO

p_pelagic

16

4

1

0.367

0.362

0.120

0.660

1.048

0.143

0.035

SSEO

p_pelagic

16

4

2

0.351

0.344

0.103

0.648

1.050

0.143

0.049

SSEO

p_pelagic

16

42

1

0.832

0.832

0.826

0.839

1.113

0.003

0.000

SSEO

p_yellow

16

42

1

0.404

0.404

0.384

0.427

1.050

0.011

0.000

SSEO

p_dsr

16

42

1

0.469

0.479

0.412

0.514

3.649

0.030

-0.330

SSEO

p_slope

16

42

1

0.145

0.132

0.106

0.204

4.992

0.032

0.411

SSEO

p_pelagic

16

5

1

0.367

0.360

0.118

0.662

1.052

0.145

0.048

SSEO

p_pelagic

16

5

2

0.350

0.343

0.112

0.641

1.048

0.140

0.050

SSEO

p_pelagic

16

6

1

0.369

0.361

0.127

0.654

1.066

0.142

0.056

SSEO

p_pelagic

16

6

2

0.348

0.342

0.107

0.631

1.062

0.140

0.043

SSEO

p_pelagic

16

7

1

0.372

0.366

0.121

0.665

1.053

0.143

0.042

SSEO

p_pelagic

16

7

2

0.349

0.344

0.109

0.634

1.056

0.140

0.036

SSEO

p_pelagic

16

8

1

0.371

0.363

0.122

0.658

1.060

0.143

0.056

SSEO

p_pelagic

16

8

2

0.346

0.338

0.113

0.642

1.042

0.140

0.057

SSEO

p_pelagic

16

9

1

0.372

0.365

0.128

0.663

1.050

0.141

0.050

SSEO

p_pelagic

16

9

2

0.355

0.349

0.115

0.647

1.048

0.141

0.042

SSEO

p_black

8

40

1

0.475

0.416

0.340

0.962

1.024

0.176

0.335

afognak

p_black

8

48

1

0.534

0.494

0.414

0.967

1.017

0.141

0.283

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

tau_beta1_pH

1

1

6.567

1.008

0.060

55.310

1.158

14.540

0.382

CI

beta3_pH

10

1

33.008

32.518

30.371

43.355

1.311

2.484

0.197

northeast

beta2_yellow

13

-7.169

-6.655

-19.459

-0.397

1.113

5.377

-0.096

NSEI

beta3_pH

16

4

41.524

41.514

40.904

42.322

1.203

0.376

0.027

SSEO

beta0_pelagic

16

-0.618

-0.495

-1.664

0.024

1.171

0.474

-0.259

SSEO

beta1_pelagic

16

1.816

1.713

0.831

3.050

1.224

0.657

0.157

SSEO

beta0_pH

2

1

1.562

1.508

1.090

2.047

2.901

0.273

0.198

NG

beta1_pH

2

1

1.329

1.443

0.432

2.399

2.891

0.561

-0.203

NG

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

418.00

183.00

1,001.00

343.00

Prop. of total

0.70

0.30

0.74

0.26

Prop. median less than 1

0.32

0.27

0.34

0.41

Prop. median less than 10

0.46

0.60

0.55

0.77

Prop. median less than 50

0.61

0.85

0.71

0.83

Median of the median values

16.59

6.73

5.84

1.80

Mean of the median values

408.74

94.17

218.02

89.32

Mean skew of posteriors

0.26

0.28

0.27

0.24

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

6.00

10.00

65.00

34.00

Prop. of total

0.38

0.62

0.66

0.34

Prop. median less than 1

1.00

0.40

0.69

0.65

Prop. median less than 10

1.00

0.90

0.82

0.94

Prop. median less than 50

1.00

1.00

0.85

0.94

Median of the median values

0.00

1.37

0.04

0.26

Mean of the median values

0.01

3.71

142.91

5.54

Mean skew of posteriors

0.17

0.21

0.16

0.17

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_pelagic

1

beta0_pH

1

beta1_pelagic

1

beta1_pH

1

beta2_yellow

1

beta3_pH

2

H_ayg

1

Hb_ay

1

Hb_ayg

1

Hb_ayu

1

Hd_ay

1

Hd_ayg

2

Ho_ay

2

Ho_ayg

1

Ho_ayu

7

Hp_ay

1

Hp_ayg

1

Hp_ayu

1

Hs_ay

2

Hs_ayg

2

Hs_ayu

3

Hy_ay

3

Hy_ayu

3

p_dsr

2

p_pelagic

1

p_slope

3

pH

19

R_ay

1

R_ayg

2

R_ayu

1

Rb_ayg

2

Rd_ay

6

Rd_ayg

6

Rd_ayu

6

re_dsr

2

re_pH

60

re_slope

2

Ro_ay

20

Ro_ayg

11

Ro_ayu

47

Rp_ayg

1

Rs_ay

9

Rs_ayg

9

Rs_ayu

19

Ry_ay

11

Ry_ayg

4

Ry_ayu

21

tau_beta1_pH

1

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.017 0.067 -0.143 -0.019 0.123
mu_bc_H[2] 0.068 0.053 -0.032 0.066 0.180
mu_bc_H[3] -0.275 0.066 -0.403 -0.278 -0.141
mu_bc_H[4] -0.931 0.182 -1.293 -0.928 -0.583
mu_bc_H[5] 1.388 1.232 0.094 1.091 4.525
mu_bc_H[6] -2.395 0.306 -2.973 -2.399 -1.780
mu_bc_H[7] -0.451 0.099 -0.656 -0.450 -0.258
mu_bc_H[8] 0.237 0.343 -0.341 0.210 1.020
mu_bc_H[9] -0.249 0.124 -0.485 -0.251 0.003
mu_bc_H[10] 0.055 0.065 -0.071 0.053 0.188
mu_bc_H[11] -0.126 0.037 -0.199 -0.125 -0.052
mu_bc_H[12] -0.249 0.107 -0.470 -0.245 -0.041
mu_bc_H[13] -0.160 0.073 -0.305 -0.162 -0.016
mu_bc_H[14] -0.308 0.090 -0.488 -0.306 -0.137
mu_bc_H[15] -0.333 0.046 -0.418 -0.334 -0.239
mu_bc_H[16] -0.156 0.371 -0.801 -0.183 0.661
mu_bc_R[1] 1.044 0.120 0.802 1.045 1.283
mu_bc_R[2] 1.018 0.075 0.867 1.021 1.162
mu_bc_R[3] 0.243 0.093 0.052 0.247 0.417
mu_bc_R[4] 0.677 0.200 0.271 0.685 1.049
mu_bc_R[5] 1.689 0.381 0.946 1.698 2.428
mu_bc_R[6] -1.436 0.334 -2.075 -1.435 -0.779
mu_bc_R[7] 0.513 0.169 0.162 0.518 0.836
mu_bc_R[8] 0.648 0.149 0.369 0.642 0.943
mu_bc_R[9] 0.401 0.167 0.048 0.409 0.704
mu_bc_R[10] 0.904 0.105 0.701 0.902 1.108
mu_bc_R[11] 1.105 0.074 0.961 1.107 1.256
mu_bc_R[12] 0.880 0.162 0.565 0.883 1.202
mu_bc_R[13] 1.042 0.091 0.859 1.045 1.218
mu_bc_R[14] 0.866 0.123 0.614 0.867 1.104
mu_bc_R[15] 0.806 0.095 0.622 0.808 0.990
mu_bc_R[16] 1.360 0.115 1.124 1.362 1.583
tau_pH[1] 2.398 0.191 2.031 2.394 2.784
tau_pH[2] 1.630 0.155 1.350 1.622 1.952
tau_pH[3] 1.217 0.113 1.008 1.212 1.454
tau_pH[4] 1.090 0.147 0.832 1.078 1.411
tau_pH[5] 0.439 0.078 0.304 0.433 0.603
beta0_pH[1,1] 0.549 0.181 0.168 0.557 0.871
beta0_pH[2,1] 1.562 0.273 1.090 1.508 2.047
beta0_pH[3,1] 1.421 0.213 0.907 1.440 1.783
beta0_pH[4,1] 1.605 0.235 1.101 1.620 2.023
beta0_pH[5,1] -0.771 0.406 -1.722 -0.723 -0.140
beta0_pH[6,1] -0.524 0.494 -1.942 -0.381 0.080
beta0_pH[7,1] 0.758 0.011 0.737 0.758 0.779
beta0_pH[8,1] -0.526 0.338 -1.308 -0.499 -0.006
beta0_pH[9,1] -0.293 0.288 -0.911 -0.272 0.211
beta0_pH[10,1] 0.250 0.241 -0.298 0.265 0.668
beta0_pH[11,1] 0.122 0.275 -0.433 0.122 0.628
beta0_pH[12,1] 0.620 0.193 0.208 0.625 0.983
beta0_pH[13,1] -0.665 0.145 -0.972 -0.659 -0.397
beta0_pH[14,1] 0.023 0.181 -0.344 0.024 0.358
beta0_pH[15,1] 0.066 0.201 -0.402 0.082 0.415
beta0_pH[16,1] -0.863 0.215 -1.319 -0.856 -0.458
beta0_pH[1,2] 2.563 0.082 2.404 2.562 2.732
beta0_pH[2,2] 3.079 0.022 3.038 3.079 3.122
beta0_pH[3,2] 3.624 0.040 3.546 3.624 3.703
beta0_pH[4,2] 2.936 0.019 2.899 2.935 2.973
beta0_pH[5,2] 4.551 1.515 2.782 4.210 8.468
beta0_pH[6,2] 3.324 0.156 3.029 3.322 3.642
beta0_pH[7,2] 2.034 0.050 1.938 2.034 2.133
beta0_pH[8,2] 2.822 0.066 2.698 2.820 2.957
beta0_pH[9,2] 3.951 0.213 3.544 3.946 4.394
beta0_pH[10,2] 3.994 0.130 3.749 3.991 4.252
beta0_pH[11,2] 1.885 0.009 1.868 1.886 1.903
beta0_pH[12,2] 1.792 0.066 1.665 1.791 1.925
beta0_pH[13,2] 2.497 0.027 2.446 2.498 2.551
beta0_pH[14,2] 1.803 0.024 1.757 1.803 1.850
beta0_pH[15,2] 2.545 0.017 2.512 2.544 2.577
beta0_pH[16,2] 2.766 0.025 2.716 2.765 2.815
beta0_pH[1,3] 1.121 0.351 0.229 1.181 1.666
beta0_pH[2,3] 2.127 0.017 2.092 2.126 2.160
beta0_pH[3,3] 2.692 0.024 2.644 2.692 2.739
beta0_pH[4,3] 2.979 0.034 2.911 2.979 3.047
beta0_pH[5,3] 3.346 2.136 0.909 2.858 8.998
beta0_pH[6,3] 1.356 0.374 0.723 1.328 2.164
beta0_pH[7,3] 0.860 0.049 0.762 0.859 0.957
beta0_pH[8,3] 0.250 0.301 -0.300 0.282 0.620
beta0_pH[9,3] -2.098 1.043 -4.509 -1.992 -0.368
beta0_pH[10,3] -0.435 0.837 -2.455 -0.327 0.877
beta0_pH[11,3] 1.096 0.511 0.027 1.130 2.026
beta0_pH[12,3] 0.947 0.652 -0.295 0.946 2.207
beta0_pH[13,3] 0.999 0.481 -0.063 1.013 1.912
beta0_pH[14,3] 0.565 0.430 -0.355 0.595 1.324
beta0_pH[15,3] 0.719 0.413 -0.180 0.742 1.487
beta0_pH[16,3] 1.185 0.606 0.014 1.161 2.404
beta0_pH[11,4] 0.863 0.461 -0.129 0.893 1.709
beta0_pH[12,4] 0.843 0.668 -0.491 0.825 2.175
beta0_pH[13,4] 0.624 0.539 -0.635 0.689 1.492
beta0_pH[14,4] 0.385 0.429 -0.536 0.415 1.134
beta0_pH[15,4] 0.535 0.527 -0.606 0.608 1.413
beta0_pH[16,4] 0.797 0.553 -0.317 0.791 1.901
beta0_pH[11,5] 2.970 0.153 2.691 2.963 3.284
beta0_pH[12,5] -0.895 0.507 -1.958 -0.883 0.037
beta0_pH[13,5] 4.896 0.413 4.116 4.876 5.767
beta0_pH[14,5] 3.878 0.363 3.204 3.854 4.633
beta0_pH[15,5] 2.288 0.122 2.054 2.284 2.533
beta0_pH[16,5] 0.269 0.068 0.168 0.259 0.435
beta1_pH[1,1] 2.569 0.469 1.817 2.510 3.775
beta1_pH[2,1] 1.329 0.561 0.432 1.443 2.399
beta1_pH[3,1] 1.755 0.425 1.070 1.709 2.832
beta1_pH[4,1] 1.834 0.414 1.123 1.798 2.750
beta1_pH[5,1] 2.299 0.529 1.467 2.244 3.596
beta1_pH[6,1] 2.085 0.791 1.048 1.883 4.169
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.349 0.917 1.179 2.100 4.559
beta1_pH[9,1] 2.049 0.676 1.161 1.920 4.020
beta1_pH[10,1] 2.161 0.514 1.332 2.111 3.398
beta1_pH[11,1] 3.268 0.460 2.488 3.215 4.344
beta1_pH[12,1] 2.910 0.286 2.365 2.910 3.474
beta1_pH[13,1] 3.474 0.245 3.020 3.466 3.975
beta1_pH[14,1] 2.814 0.324 2.209 2.802 3.504
beta1_pH[15,1] 2.383 0.290 1.870 2.363 3.045
beta1_pH[16,1] 4.015 0.299 3.450 4.016 4.624
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 1.904 1.116 0.076 1.783 4.706
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.714 2.406 0.792 4.256 10.812
beta1_pH[9,3] 7.257 2.671 3.242 6.817 13.412
beta1_pH[10,3] 6.326 2.587 2.736 5.768 12.409
beta1_pH[11,3] 2.198 0.552 1.162 2.194 3.311
beta1_pH[12,3] 3.274 0.663 2.049 3.256 4.625
beta1_pH[13,3] 1.999 0.513 0.988 1.997 3.076
beta1_pH[14,3] 2.219 0.485 1.333 2.194 3.218
beta1_pH[15,3] 2.107 0.478 1.220 2.084 3.115
beta1_pH[16,3] 3.703 0.660 2.369 3.730 4.982
beta1_pH[11,4] 1.940 0.514 0.941 1.932 2.959
beta1_pH[12,4] 3.210 0.676 1.885 3.199 4.617
beta1_pH[13,4] 1.384 0.554 0.427 1.336 2.645
beta1_pH[14,4] 1.994 0.497 1.068 1.975 3.030
beta1_pH[15,4] 1.472 0.568 0.477 1.435 2.644
beta1_pH[16,4] 3.670 0.641 2.366 3.691 4.935
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.244 0.838 0.335 0.942 3.407
beta2_pH[2,1] 1.780 0.970 0.297 1.686 3.448
beta2_pH[3,1] 1.445 0.838 0.224 1.325 2.959
beta2_pH[4,1] 1.229 0.828 0.242 0.951 2.954
beta2_pH[5,1] 2.322 2.769 0.264 1.466 10.454
beta2_pH[6,1] 1.504 2.765 0.114 0.476 9.699
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 1.314 2.367 0.128 0.417 9.473
beta2_pH[9,1] 1.485 2.427 0.143 0.580 9.662
beta2_pH[10,1] 1.413 2.417 -1.725 0.778 9.486
beta2_pH[11,1] 0.437 0.239 0.195 0.376 0.999
beta2_pH[12,1] 1.550 0.811 0.595 1.286 3.624
beta2_pH[13,1] 1.031 0.454 0.543 0.917 2.309
beta2_pH[14,1] 1.170 0.689 0.449 0.957 3.272
beta2_pH[15,1] 1.061 0.654 0.353 0.873 3.044
beta2_pH[16,1] 1.055 0.331 0.628 0.993 1.858
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 7.024 6.875 0.167 5.548 22.760
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 1.391 1.099 0.129 1.066 4.203
beta2_pH[9,3] 0.307 0.492 0.089 0.202 1.187
beta2_pH[10,3] 0.442 0.690 0.109 0.270 2.006
beta2_pH[11,3] -5.869 4.595 -17.625 -4.626 -0.592
beta2_pH[12,3] -6.640 4.498 -17.968 -5.585 -1.229
beta2_pH[13,3] -6.090 4.535 -18.011 -4.758 -1.074
beta2_pH[14,3] -6.066 4.436 -18.014 -4.771 -1.163
beta2_pH[15,3] -6.263 4.347 -17.555 -5.079 -1.116
beta2_pH[16,3] -6.257 4.375 -17.881 -4.919 -1.546
beta2_pH[11,4] -6.192 4.464 -17.347 -5.036 -0.797
beta2_pH[12,4] -6.493 4.502 -17.701 -5.301 -1.161
beta2_pH[13,4] -6.184 4.459 -17.266 -4.999 -0.847
beta2_pH[14,4] -6.280 4.370 -17.660 -5.056 -1.215
beta2_pH[15,4] -6.246 4.529 -17.669 -5.021 -0.840
beta2_pH[16,4] -6.654 4.217 -17.378 -5.581 -1.648
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.922 1.122 32.834 34.928 37.362
beta3_pH[2,1] 32.715 1.597 29.633 32.660 36.202
beta3_pH[3,1] 33.006 1.517 30.083 32.979 36.497
beta3_pH[4,1] 32.640 1.779 29.512 32.552 36.365
beta3_pH[5,1] 27.692 1.225 25.970 27.531 31.065
beta3_pH[6,1] 32.978 3.145 25.607 33.112 39.792
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 35.885 3.851 28.932 36.513 42.231
beta3_pH[9,1] 32.490 2.967 27.755 32.000 39.785
beta3_pH[10,1] 33.008 2.484 30.371 32.518 43.355
beta3_pH[11,1] 31.393 1.144 29.366 31.375 33.636
beta3_pH[12,1] 30.780 0.492 29.776 30.805 31.744
beta3_pH[13,1] 31.423 0.548 30.363 31.408 32.581
beta3_pH[14,1] 32.471 0.738 31.233 32.416 34.068
beta3_pH[15,1] 30.974 0.726 29.467 31.006 32.321
beta3_pH[16,1] 30.127 0.354 29.423 30.132 30.818
beta3_pH[1,2] 34.428 5.546 25.500 34.477 43.545
beta3_pH[2,2] 34.479 5.521 25.438 34.607 43.471
beta3_pH[3,2] 30.258 7.239 18.619 30.064 42.284
beta3_pH[4,2] 30.501 7.101 18.681 30.557 42.277
beta3_pH[5,2] 30.434 7.217 18.563 30.349 42.406
beta3_pH[6,2] 30.368 7.259 18.619 30.376 42.212
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 38.446 3.219 30.512 39.648 42.669
beta3_pH[2,3] 34.530 5.466 25.432 34.617 43.515
beta3_pH[3,3] 30.430 7.063 18.713 30.345 42.427
beta3_pH[4,3] 30.589 7.195 18.619 31.004 42.292
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.023 1.321 38.503 42.073 43.857
beta3_pH[9,3] 36.203 1.944 34.521 35.423 41.856
beta3_pH[10,3] 37.788 2.502 34.632 37.083 43.286
beta3_pH[11,3] 40.722 1.077 38.909 40.735 42.719
beta3_pH[12,3] 42.019 0.730 40.310 42.258 42.851
beta3_pH[13,3] 41.361 0.720 39.931 41.364 42.730
beta3_pH[14,3] 40.488 0.669 39.130 40.523 41.780
beta3_pH[15,3] 40.637 0.642 39.160 40.654 41.849
beta3_pH[16,3] 41.840 0.373 41.030 41.903 42.548
beta3_pH[11,4] 40.557 1.079 38.719 40.451 42.993
beta3_pH[12,4] 41.898 0.791 40.288 42.150 42.908
beta3_pH[13,4] 41.868 1.102 39.721 41.929 43.780
beta3_pH[14,4] 40.434 0.712 38.973 40.467 41.725
beta3_pH[15,4] 41.317 1.085 39.045 41.269 43.384
beta3_pH[16,4] 41.524 0.376 40.904 41.514 42.322
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.264 0.021 0.224 0.264 0.305
beta4_pH[2,1] 0.211 0.005 0.201 0.211 0.221
beta4_pH[3,1] 0.172 0.007 0.159 0.172 0.186
beta4_pH[4,1] 0.294 0.011 0.271 0.294 0.316
beta4_pH[5,1] 0.422 0.271 0.092 0.388 0.937
beta4_pH[6,1] 0.787 0.382 0.343 0.738 1.484
beta4_pH[7,1] 0.843 0.335 0.351 0.794 1.609
beta4_pH[8,1] 0.161 0.071 0.046 0.155 0.313
beta4_pH[9,1] 0.217 0.087 0.062 0.213 0.400
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.447 0.011 0.425 0.447 0.469
beta4_pH[12,1] 0.642 0.023 0.599 0.642 0.687
beta4_pH[13,1] 0.630 0.014 0.604 0.630 0.658
beta4_pH[14,1] 0.671 0.032 0.608 0.671 0.734
beta4_pH[15,1] 0.799 0.021 0.759 0.799 0.840
beta4_pH[16,1] 1.224 0.043 1.142 1.222 1.314
beta4_pH[1,2] 0.306 0.097 0.133 0.301 0.514
beta4_pH[2,2] 0.514 0.023 0.470 0.514 0.561
beta4_pH[3,2] 0.446 0.012 0.423 0.446 0.469
beta4_pH[4,2] 0.654 0.037 0.582 0.654 0.728
beta4_pH[5,2] 0.747 0.950 0.111 0.536 2.603
beta4_pH[6,2] 0.911 0.731 0.219 0.742 2.654
beta4_pH[7,2] 0.880 0.404 0.344 0.803 1.922
beta4_pH[8,2] 0.846 0.205 0.482 0.832 1.270
beta4_pH[9,2] 1.152 1.057 0.304 0.883 3.566
beta4_pH[10,2] 0.454 0.051 0.359 0.454 0.558
beta4_pH[11,2] 1.269 0.050 1.174 1.268 1.369
beta4_pH[12,2] 0.580 0.108 0.382 0.574 0.802
beta4_pH[13,2] 0.945 0.045 0.859 0.945 1.035
beta4_pH[14,2] 0.931 0.061 0.812 0.930 1.055
beta4_pH[15,2] 1.010 0.026 0.960 1.010 1.061
beta4_pH[16,2] 1.407 0.074 1.264 1.407 1.555
beta4_pH[1,3] 0.283 0.088 0.123 0.280 0.470
beta4_pH[2,3] 0.051 0.017 0.021 0.050 0.086
beta4_pH[3,3] 0.052 0.010 0.034 0.052 0.071
beta4_pH[4,3] 0.174 0.032 0.115 0.173 0.238
beta4_pH[5,3] 0.730 0.738 0.113 0.566 2.276
beta4_pH[6,3] 0.856 0.659 0.209 0.714 2.358
beta4_pH[7,3] 0.886 0.403 0.344 0.817 1.854
beta4_pH[8,3] 0.406 0.189 0.107 0.387 0.842
beta4_pH[9,3] 0.949 1.065 0.176 0.720 2.950
beta4_pH[10,3] 0.061 0.041 0.003 0.055 0.154
beta4_pH[11,3] 0.364 0.016 0.332 0.363 0.394
beta4_pH[12,3] 0.440 0.042 0.360 0.439 0.526
beta4_pH[13,3] 0.836 0.037 0.767 0.835 0.909
beta4_pH[14,3] 0.824 0.067 0.700 0.822 0.962
beta4_pH[15,3] 0.635 0.017 0.603 0.634 0.668
beta4_pH[16,3] 0.635 0.035 0.570 0.634 0.706
beta4_pH[1,4] 0.644 0.347 0.202 0.577 1.519
beta4_pH[2,4] 0.807 0.844 0.102 0.586 2.838
beta4_pH[3,4] 0.859 0.904 0.117 0.609 3.081
beta4_pH[4,4] 0.749 0.494 0.215 0.635 1.932
beta4_pH[5,4] 0.770 0.857 0.108 0.566 2.652
beta4_pH[6,4] 0.917 0.866 0.209 0.721 2.734
beta4_pH[7,4] 0.891 0.425 0.348 0.816 1.858
beta4_pH[8,4] 0.722 0.601 0.152 0.579 2.152
beta4_pH[9,4] 0.791 0.932 0.128 0.565 2.846
beta4_pH[10,4] 1.324 2.753 0.007 0.372 9.666
beta4_pH[11,4] 0.444 0.017 0.410 0.444 0.478
beta4_pH[12,4] 0.557 0.054 0.455 0.556 0.663
beta4_pH[13,4] 0.925 0.036 0.856 0.924 0.999
beta4_pH[14,4] 0.919 0.073 0.782 0.918 1.062
beta4_pH[15,4] 0.714 0.020 0.675 0.714 0.754
beta4_pH[16,4] 0.684 0.041 0.607 0.683 0.769
beta4_pH[1,5] 0.624 0.331 0.215 0.556 1.421
beta4_pH[2,5] 0.782 0.831 0.102 0.560 2.734
beta4_pH[3,5] 0.857 0.982 0.110 0.600 3.210
beta4_pH[4,5] 0.717 0.437 0.199 0.616 1.793
beta4_pH[5,5] 0.775 0.833 0.114 0.587 2.666
beta4_pH[6,5] 0.906 0.759 0.215 0.734 2.756
beta4_pH[7,5] 0.880 0.436 0.338 0.805 1.923
beta4_pH[8,5] 0.742 0.634 0.157 0.591 2.253
beta4_pH[9,5] 0.778 0.890 0.111 0.562 2.682
beta4_pH[10,5] 1.477 2.967 0.007 0.395 10.326
beta4_pH[11,5] 0.604 0.051 0.507 0.604 0.706
beta4_pH[12,5] 0.585 0.167 0.308 0.567 0.966
beta4_pH[13,5] 0.774 0.062 0.661 0.771 0.905
beta4_pH[14,5] 0.843 0.116 0.634 0.839 1.087
beta4_pH[15,5] 1.141 0.072 1.006 1.140 1.286
beta4_pH[16,5] 17.693 4.055 9.966 17.689 24.484
beta0_pelagic[1] 2.268 0.036 2.200 2.267 2.340
beta0_pelagic[2] 1.537 0.007 1.524 1.538 1.551
beta0_pelagic[3] 0.196 0.234 -0.338 0.243 0.560
beta0_pelagic[4] 0.230 0.218 -0.311 0.248 0.601
beta0_pelagic[5] 1.446 0.155 1.148 1.443 1.749
beta0_pelagic[6] 1.354 0.060 1.237 1.354 1.476
beta0_pelagic[7] 1.483 0.038 1.414 1.481 1.562
beta0_pelagic[8] 1.675 0.043 1.596 1.673 1.764
beta0_pelagic[9] 2.793 0.065 2.668 2.791 2.927
beta0_pelagic[10] 2.602 0.032 2.543 2.602 2.667
beta0_pelagic[11] -0.136 0.164 -0.501 -0.121 0.153
beta0_pelagic[12] 1.801 0.036 1.735 1.800 1.874
beta0_pelagic[13] -0.835 0.449 -1.661 -0.836 -0.004
beta0_pelagic[14] -0.280 0.176 -0.667 -0.264 0.018
beta0_pelagic[15] -0.215 0.039 -0.317 -0.198 -0.181
beta0_pelagic[16] -0.618 0.474 -1.664 -0.495 0.024
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.749 0.251 0.318 0.728 1.345
beta1_pelagic[4] 0.818 0.243 0.395 0.802 1.408
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.150 0.235 0.741 1.130 1.698
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 2.574 0.578 1.323 2.635 3.451
beta1_pelagic[14] 1.784 0.443 1.107 1.712 2.921
beta1_pelagic[15] 0.130 0.197 0.000 0.021 0.669
beta1_pelagic[16] 1.816 0.657 0.831 1.713 3.050
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 6.378 5.958 0.233 4.630 21.151
beta2_pelagic[4] 7.503 5.921 0.356 6.311 21.346
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.807 1.527 0.276 1.291 6.103
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.201 0.102 0.104 0.175 0.472
beta2_pelagic[14] 0.884 0.945 0.202 0.585 3.764
beta2_pelagic[15] 2.632 1.870 0.170 2.331 7.072
beta2_pelagic[16] 0.603 0.913 0.111 0.285 3.734
beta3_pelagic[1] 28.733 3.855 22.406 28.660 35.187
beta3_pelagic[2] 28.745 3.905 22.276 28.746 35.197
beta3_pelagic[3] 27.789 2.705 23.306 28.031 32.161
beta3_pelagic[4] 26.113 1.303 23.330 26.272 28.833
beta3_pelagic[5] 29.494 6.707 18.659 29.563 40.311
beta3_pelagic[6] 29.596 6.596 18.559 29.715 40.409
beta3_pelagic[7] 29.588 6.606 18.655 29.763 40.382
beta3_pelagic[8] 29.743 6.524 18.591 29.999 40.408
beta3_pelagic[9] 29.461 6.666 18.597 29.330 40.377
beta3_pelagic[10] 29.638 6.580 18.678 29.782 40.448
beta3_pelagic[11] 30.204 1.116 28.121 30.158 32.513
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 31.661 3.873 25.181 31.370 39.688
beta3_pelagic[14] 34.243 1.710 31.522 34.027 38.325
beta3_pelagic[15] 34.636 5.778 23.731 35.949 41.749
beta3_pelagic[16] 30.891 3.268 24.789 30.901 38.281
beta0_yellow[1] -0.482 0.078 -0.641 -0.484 -0.326
beta0_yellow[2] 0.547 0.080 0.360 0.551 0.687
beta0_yellow[3] -0.214 0.129 -0.480 -0.219 0.013
beta0_yellow[4] 0.887 0.248 0.222 0.948 1.192
beta0_yellow[5] 0.260 0.322 -0.424 0.283 0.840
beta0_yellow[6] 0.819 0.106 0.614 0.819 1.029
beta0_yellow[7] 0.996 0.095 0.805 0.994 1.179
beta0_yellow[8] 1.026 0.094 0.840 1.026 1.205
beta0_yellow[9] 0.459 0.103 0.260 0.457 0.658
beta0_yellow[10] 0.537 0.070 0.400 0.537 0.675
beta0_yellow[11] -0.192 0.125 -0.457 -0.182 0.019
beta0_yellow[12] -1.082 0.093 -1.265 -1.081 -0.902
beta0_yellow[13] -1.348 0.160 -1.689 -1.337 -1.069
beta0_yellow[14] -0.253 0.172 -0.580 -0.233 -0.014
beta0_yellow[15] -1.256 0.147 -1.623 -1.232 -1.038
beta0_yellow[16] -0.252 0.092 -0.430 -0.254 -0.069
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.006 0.237 0.559 1.002 1.476
beta1_yellow[3] 0.478 0.240 0.022 0.485 0.921
beta1_yellow[4] 1.339 0.594 0.634 1.190 3.202
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.307 0.404 0.001 0.093 1.293
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 1.089 1.188 0.025 0.757 5.018
beta1_yellow[14] 0.742 1.069 0.013 0.415 4.327
beta1_yellow[15] 0.166 0.197 0.011 0.082 0.689
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -7.795 5.781 -21.678 -6.575 -0.481
beta2_yellow[3] -7.590 5.686 -21.402 -6.383 -0.396
beta2_yellow[4] -5.500 5.808 -20.065 -3.599 -0.128
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -7.969 5.913 -22.590 -6.759 -0.429
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -7.169 5.377 -19.459 -6.655 -0.397
beta2_yellow[14] -7.320 5.555 -20.074 -6.003 -0.381
beta2_yellow[15] -7.719 5.722 -21.233 -6.573 -0.364
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.852 1.110 27.329 28.651 31.764
beta3_yellow[3] 32.410 3.231 23.338 32.572 38.398
beta3_yellow[4] 29.946 3.300 23.867 28.719 35.769
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 34.365 3.523 29.572 33.470 42.360
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 33.795 2.804 29.222 33.908 40.483
beta3_yellow[14] 33.955 3.643 29.152 33.529 42.971
beta3_yellow[15] 36.649 3.881 30.343 36.458 43.587
beta3_yellow[16] 25.000 0.000 25.000 25.000 25.000
beta0_black[1] 0.011 0.075 -0.138 0.012 0.157
beta0_black[2] 1.855 0.031 1.795 1.855 1.919
beta0_black[3] 1.238 0.037 1.168 1.238 1.314
beta0_black[4] 2.424 0.048 2.333 2.422 2.520
beta0_black[5] 5.130 2.408 1.843 4.606 11.141
beta0_black[6] 5.212 2.225 2.626 4.626 11.187
beta0_black[7] 4.150 2.181 1.583 3.565 10.091
beta0_black[8] 0.879 0.215 0.483 0.870 1.328
beta0_black[9] 2.798 0.269 2.277 2.799 3.322
beta0_black[10] 1.287 0.052 1.183 1.289 1.386
beta0_black[11] 2.978 0.052 2.879 2.978 3.081
beta0_black[12] 4.913 0.192 4.532 4.910 5.297
beta0_black[13] 0.429 0.329 0.018 0.353 1.223
beta0_black[14] 2.853 0.106 2.637 2.855 3.054
beta0_black[15] 0.723 0.051 0.625 0.723 0.823
beta0_black[16] 3.772 0.072 3.629 3.772 3.913
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.545 1.803 -6.833 -2.136 -0.295
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.440 3.136 28.044 38.165 40.884
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.325 0.092 -0.500 -0.326 -0.139
beta4_black[2] 0.258 0.052 0.154 0.259 0.362
beta4_black[3] -0.719 0.064 -0.848 -0.720 -0.595
beta4_black[4] 0.632 0.127 0.386 0.634 0.884
beta4_black[5] 0.704 1.302 -1.079 0.460 3.843
beta4_black[6] 0.847 1.401 -0.826 0.559 4.261
beta4_black[7] 0.727 1.295 -0.953 0.483 3.969
beta4_black[8] 0.089 0.276 -0.476 0.087 0.628
beta4_black[9] 1.028 0.830 -0.080 0.859 3.052
beta4_black[10] 0.073 0.067 -0.056 0.075 0.203
beta4_black[11] -0.608 0.064 -0.731 -0.608 -0.485
beta4_black[12] 1.234 0.294 0.686 1.225 1.853
beta4_black[13] -1.004 0.077 -1.152 -1.003 -0.856
beta4_black[14] 0.220 0.184 -0.140 0.218 0.592
beta4_black[15] 0.002 0.046 -0.089 0.002 0.093
beta4_black[16] -0.416 0.117 -0.647 -0.417 -0.184
beta0_dsr[11] -0.329 0.010 -0.348 -0.328 -0.310
beta0_dsr[12] 1.409 0.043 1.327 1.408 1.496
beta0_dsr[13] 0.155 0.022 0.112 0.154 0.198
beta0_dsr[14] -0.126 0.017 -0.160 -0.127 -0.093
beta0_dsr[15] 0.206 0.014 0.178 0.206 0.233
beta0_dsr[16] 0.166 0.016 0.134 0.166 0.197
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.544 0.036 0.475 0.544 0.613
beta4_dsr[12] 0.799 0.226 0.366 0.795 1.249
beta4_dsr[13] -0.078 0.041 -0.160 -0.078 0.002
beta4_dsr[14] -0.208 0.071 -0.350 -0.207 -0.071
beta4_dsr[15] 0.216 0.022 0.174 0.216 0.259
beta4_dsr[16] -0.055 0.045 -0.142 -0.055 0.033
beta0_slope[11] -2.209 0.021 -2.249 -2.209 -2.168
beta0_slope[12] -4.984 0.192 -5.381 -4.984 -4.622
beta0_slope[13] -1.500 0.034 -1.568 -1.500 -1.433
beta0_slope[14] -2.712 0.051 -2.813 -2.710 -2.616
beta0_slope[15] -1.592 0.022 -1.637 -1.593 -1.549
beta0_slope[16] -2.661 0.037 -2.730 -2.661 -2.589
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.333 0.076 -0.489 -0.332 -0.183
beta4_slope[12] -4.658 0.815 -6.395 -4.614 -3.124
beta4_slope[13] 0.326 0.053 0.221 0.326 0.434
beta4_slope[14] -0.161 0.152 -0.472 -0.158 0.133
beta4_slope[15] -0.360 0.034 -0.429 -0.360 -0.296
beta4_slope[16] -0.136 0.102 -0.339 -0.133 0.063
beta5_slope[11] 3.065 0.104 2.870 3.061 3.270
beta5_slope[12] 2.005 0.390 1.276 1.989 2.829
beta5_slope[13] 0.877 0.078 0.726 0.880 1.029
beta5_slope[14] 3.754 0.352 3.107 3.736 4.496
beta5_slope[15] 1.305 0.072 1.166 1.304 1.445
beta5_slope[16] 2.728 0.237 2.291 2.720 3.220
tau_comp[1] 3.801 0.302 3.260 3.795 4.429
tau_comp[2] 2.863 0.256 2.404 2.847 3.396
tau_comp[3] 0.871 0.073 0.730 0.871 1.019
tau_comp[4] 1.793 0.274 1.314 1.777 2.378
tau_comp[5] 1.290 0.187 0.958 1.282 1.670
sigma_H[1] 0.208 0.049 0.121 0.204 0.310
sigma_H[2] 0.175 0.029 0.127 0.172 0.238
sigma_H[3] 0.187 0.040 0.114 0.185 0.269
sigma_H[4] 0.458 0.080 0.324 0.449 0.640
sigma_H[5] 0.960 0.203 0.599 0.948 1.395
sigma_H[6] 0.534 0.191 0.149 0.529 0.914
sigma_H[7] 0.310 0.061 0.215 0.303 0.454
sigma_H[8] 0.419 0.099 0.253 0.411 0.622
sigma_H[9] 0.510 0.116 0.327 0.494 0.778
sigma_H[10] 0.211 0.039 0.144 0.209 0.296
sigma_H[11] 0.266 0.043 0.194 0.263 0.365
sigma_H[12] 0.379 0.145 0.198 0.336 0.713
sigma_H[13] 0.217 0.036 0.156 0.214 0.297
sigma_H[14] 0.499 0.089 0.349 0.491 0.691
sigma_H[15] 0.231 0.038 0.167 0.228 0.313
sigma_H[16] 0.214 0.041 0.149 0.209 0.309
lambda_H[1] 2.133 3.227 0.076 1.021 11.178
lambda_H[2] 5.103 4.738 0.466 3.697 17.823
lambda_H[3] 5.491 7.875 0.281 2.914 27.130
lambda_H[4] 0.007 0.004 0.001 0.005 0.018
lambda_H[5] 2.674 6.678 0.023 0.695 20.697
lambda_H[6] 5.894 12.680 0.008 0.299 43.591
lambda_H[7] 0.012 0.009 0.002 0.010 0.034
lambda_H[8] 7.608 9.920 0.096 3.921 37.733
lambda_H[9] 0.015 0.010 0.003 0.013 0.041
lambda_H[10] 0.225 0.275 0.024 0.152 0.784
lambda_H[11] 0.234 0.345 0.010 0.114 1.116
lambda_H[12] 3.826 5.236 0.206 2.211 17.805
lambda_H[13] 3.110 2.811 0.260 2.285 10.419
lambda_H[14] 3.243 3.777 0.205 2.067 13.411
lambda_H[15] 0.018 0.018 0.003 0.014 0.060
lambda_H[16] 0.559 0.981 0.031 0.271 2.883
mu_lambda_H[1] 4.094 1.900 1.017 3.867 8.243
mu_lambda_H[2] 3.708 1.940 0.504 3.566 7.822
mu_lambda_H[3] 3.375 1.826 0.750 3.093 7.640
sigma_lambda_H[1] 8.152 4.340 1.695 7.501 18.254
sigma_lambda_H[2] 8.126 4.668 0.813 7.620 18.251
sigma_lambda_H[3] 6.231 3.989 1.042 5.397 16.133
beta_H[1,1] 6.489 1.332 3.299 6.736 8.474
beta_H[2,1] 9.732 0.637 8.295 9.793 10.833
beta_H[3,1] 7.840 0.783 5.998 7.942 9.164
beta_H[4,1] 10.621 7.704 -4.807 10.616 25.554
beta_H[5,1] -0.306 2.603 -5.689 -0.183 4.171
beta_H[6,1] 3.067 4.363 -7.688 4.577 8.202
beta_H[7,1] 1.103 5.642 -11.397 1.414 11.100
beta_H[8,1] 1.320 3.978 -2.654 1.235 3.518
beta_H[9,1] 13.662 5.719 2.514 13.549 25.445
beta_H[10,1] 6.805 1.840 2.942 6.839 10.363
beta_H[11,1] 4.871 3.549 -2.983 5.549 9.899
beta_H[12,1] 2.585 1.059 0.651 2.521 4.940
beta_H[13,1] 9.072 0.915 7.147 9.137 10.625
beta_H[14,1] 2.135 1.041 0.020 2.156 4.164
beta_H[15,1] -7.046 3.344 -13.310 -7.189 -0.078
beta_H[16,1] 3.939 3.007 -0.750 3.479 10.906
beta_H[1,2] 7.773 0.269 7.228 7.776 8.274
beta_H[2,2] 9.884 0.165 9.567 9.885 10.200
beta_H[3,2] 8.801 0.198 8.416 8.801 9.209
beta_H[4,2] 3.013 1.470 0.153 2.980 6.037
beta_H[5,2] 1.757 0.998 -0.154 1.752 3.752
beta_H[6,2] 6.028 1.071 3.570 6.192 7.691
beta_H[7,2] 2.409 1.063 0.503 2.336 4.699
beta_H[8,2] 3.002 1.123 1.288 3.156 4.251
beta_H[9,2] 3.131 1.097 1.056 3.101 5.422
beta_H[10,2] 8.036 0.369 7.253 8.042 8.748
beta_H[11,2] 9.822 0.639 8.862 9.710 11.256
beta_H[12,2] 3.919 0.362 3.235 3.910 4.676
beta_H[13,2] 9.157 0.244 8.715 9.146 9.639
beta_H[14,2] 4.029 0.344 3.366 4.032 4.726
beta_H[15,2] 11.548 0.608 10.288 11.573 12.705
beta_H[16,2] 4.458 0.812 2.965 4.440 6.118
beta_H[1,3] 8.455 0.266 7.981 8.445 8.986
beta_H[2,3] 9.881 0.151 9.583 9.881 10.179
beta_H[3,3] 9.471 0.163 9.146 9.469 9.800
beta_H[4,3] -2.107 0.862 -3.796 -2.082 -0.528
beta_H[5,3] 3.890 0.620 2.579 3.901 5.082
beta_H[6,3] 8.355 1.161 6.670 8.057 10.793
beta_H[7,3] -2.585 0.626 -3.785 -2.578 -1.369
beta_H[8,3] 5.321 0.540 4.700 5.241 6.395
beta_H[9,3] -2.663 0.691 -4.107 -2.661 -1.355
beta_H[10,3] 8.571 0.278 8.028 8.570 9.134
beta_H[11,3] 8.504 0.291 7.871 8.523 9.021
beta_H[12,3] 5.284 0.304 4.565 5.320 5.792
beta_H[13,3] 8.861 0.179 8.488 8.866 9.207
beta_H[14,3] 5.786 0.268 5.188 5.807 6.245
beta_H[15,3] 10.274 0.299 9.679 10.270 10.878
beta_H[16,3] 5.963 0.614 4.736 5.987 7.092
beta_H[1,4] 8.096 0.184 7.724 8.107 8.429
beta_H[2,4] 9.861 0.130 9.590 9.868 10.103
beta_H[3,4] 9.968 0.155 9.624 9.981 10.248
beta_H[4,4] 11.799 0.461 10.893 11.807 12.680
beta_H[5,4] 5.582 0.737 4.401 5.504 7.229
beta_H[6,4] 7.194 0.901 5.310 7.441 8.458
beta_H[7,4] 8.541 0.347 7.820 8.546 9.219
beta_H[8,4] 6.735 0.269 6.200 6.741 7.205
beta_H[9,4] 7.422 0.450 6.546 7.406 8.360
beta_H[10,4] 7.530 0.226 7.112 7.525 7.985
beta_H[11,4] 9.507 0.191 9.126 9.510 9.875
beta_H[12,4] 7.241 0.192 6.860 7.242 7.626
beta_H[13,4] 9.066 0.141 8.770 9.070 9.331
beta_H[14,4] 7.834 0.217 7.416 7.833 8.279
beta_H[15,4] 9.472 0.225 9.034 9.475 9.918
beta_H[16,4] 9.414 0.247 8.960 9.411 9.892
beta_H[1,5] 8.898 0.159 8.567 8.904 9.199
beta_H[2,5] 10.672 0.101 10.479 10.669 10.879
beta_H[3,5] 10.761 0.155 10.480 10.755 11.076
beta_H[4,5] 8.352 0.486 7.410 8.345 9.327
beta_H[5,5] 5.178 0.631 3.611 5.260 6.201
beta_H[6,5] 8.969 0.636 8.031 8.839 10.361
beta_H[7,5] 6.570 0.349 5.882 6.561 7.262
beta_H[8,5] 8.251 0.214 7.884 8.238 8.698
beta_H[9,5] 8.199 0.444 7.284 8.205 9.099
beta_H[10,5] 10.039 0.227 9.579 10.044 10.492
beta_H[11,5] 11.594 0.221 11.175 11.591 12.032
beta_H[12,5] 8.604 0.182 8.263 8.603 8.972
beta_H[13,5] 10.142 0.129 9.901 10.140 10.397
beta_H[14,5] 9.249 0.218 8.844 9.243 9.692
beta_H[15,5] 11.306 0.238 10.840 11.309 11.768
beta_H[16,5] 9.836 0.181 9.479 9.841 10.180
beta_H[1,6] 10.127 0.244 9.751 10.091 10.738
beta_H[2,6] 11.473 0.114 11.251 11.469 11.702
beta_H[3,6] 10.616 0.157 10.299 10.623 10.903
beta_H[4,6] 12.906 0.872 11.178 12.926 14.571
beta_H[5,6] 5.682 0.658 4.583 5.644 7.177
beta_H[6,6] 8.740 0.580 7.442 8.786 9.834
beta_H[7,6] 10.541 0.607 9.361 10.554 11.752
beta_H[8,6] 9.524 0.274 9.038 9.541 9.972
beta_H[9,6] 8.272 0.763 6.825 8.267 9.866
beta_H[10,6] 9.142 0.346 8.411 9.169 9.758
beta_H[11,6] 10.700 0.363 9.935 10.730 11.352
beta_H[12,6] 9.371 0.224 8.948 9.358 9.845
beta_H[13,6] 11.095 0.151 10.807 11.093 11.418
beta_H[14,6] 9.915 0.282 9.341 9.925 10.447
beta_H[15,6] 10.560 0.418 9.755 10.564 11.362
beta_H[16,6] 10.592 0.250 10.054 10.603 11.076
beta_H[1,7] 10.127 1.270 6.743 10.415 11.856
beta_H[2,7] 12.113 0.508 10.995 12.153 13.025
beta_H[3,7] 10.290 0.660 8.842 10.350 11.444
beta_H[4,7] 1.426 4.449 -7.072 1.389 10.292
beta_H[5,7] 5.738 2.168 1.452 5.849 9.923
beta_H[6,7] 7.201 3.213 -0.883 8.336 11.071
beta_H[7,7] 7.669 3.063 1.421 7.663 13.700
beta_H[8,7] 10.885 1.073 9.171 10.853 12.648
beta_H[9,7] 5.025 3.979 -3.139 5.126 12.598
beta_H[10,7] 10.016 1.637 7.196 9.831 13.668
beta_H[11,7] 11.325 1.783 8.115 11.119 15.320
beta_H[12,7] 9.799 0.917 7.730 9.874 11.379
beta_H[13,7] 11.703 0.714 10.058 11.763 12.908
beta_H[14,7] 10.626 0.913 8.766 10.653 12.317
beta_H[15,7] 13.885 2.186 9.551 13.894 18.123
beta_H[16,7] 12.849 1.367 10.483 12.738 15.863
beta0_H[1] 8.213 18.251 -28.067 8.656 44.233
beta0_H[2] 10.739 8.137 -5.417 10.688 27.854
beta0_H[3] 9.823 9.884 -9.756 9.873 30.200
beta0_H[4] 3.950 177.697 -374.155 6.330 374.062
beta0_H[5] 4.550 28.430 -48.260 4.154 63.775
beta0_H[6] 5.954 52.749 -107.586 7.264 121.956
beta0_H[7] 5.853 140.816 -281.466 3.489 289.705
beta0_H[8] 7.065 37.298 -18.636 6.493 28.961
beta0_H[9] 6.805 122.824 -236.044 3.880 256.737
beta0_H[10] 8.013 38.091 -70.290 7.958 83.546
beta0_H[11] 9.654 51.823 -106.312 9.941 118.762
beta0_H[12] 6.495 11.607 -17.296 6.385 29.894
beta0_H[13] 10.167 11.895 -9.574 10.165 31.793
beta0_H[14] 7.499 11.384 -16.418 7.419 30.988
beta0_H[15] 8.150 114.970 -221.425 7.547 239.476
beta0_H[16] 7.681 29.583 -48.506 7.573 71.512